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Abstract 

I review recent developments in the CGC approach to high-energy colHsions. The focus is on 
topics related to the Quark Matter conference, specifically on predictions for the p+Pb run at the 
LHC; as an added bonus some of these predictions are confronted with first data taken during the 
p+Pb machine performance test at the LHC performed in September 2012. 1 also highlight recent 
work related to particle production fluctuations, suppression of nuclear modification factors in 
p+A from rcBK evolution, (partial) NLO corrections in forward production, and photon-hadron 
' as well as di-hadron correlations. 
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D . 1. Introduction 



In the Bjorken-Feynman parton model the constituents of a boosted hadron are viewed as 
dilute non-interacting quanta coUinear to the beam. This picture has had tremendous success 
describing the proton structure at very short distances, high-Ej jets in hadronic collisions, and 
more. However, for a probe of fixed size r (well below the confinement scale of ~ 1 fm) in the 
high-energy limit this picture is expected to break down. This is due to ample radiation in QCD 
of soft gluons as phase space opens up [1]. As a consequence, although the QCD coupling 
at the scale r may be small there are many (9(q'^') partners to interact with. Thus, at high 
energies hadrons rather correspond to a dense system of gluons and non-perturbatively strong 
fields, A'^ ~ 1/^ Ii2]- Thanks to the high occupation number these can be described as classical 
color fields sourced by the "valence" charges at higher rapidities which have been integrated 
out ifsl- The Color Glass Condensate (CGC) is an effective theory describing the dynamics of 
non-linear color fields at high energies. Exploring its properties and systematically improving its 
accuracy is one of the main goals of the high-energy nuclear physics programme. 

At present, the most suitable process for studies of high gluon density QCD is p-i-A collisions 
where a relatively dilute projectile probes a dense target. Particle production in the forward 
region of d+Au collisions at RHIC energy has discovered an interesting suppression relative to 
a simple superposition of binary p-i-p collisions fT]. However, due to limited phase space those 
measurements are restricted to relatively low transverse momenta and can not probe accurately 
the tails of the intrinsic gluon transverse momentum distributions. Also, quantitative calculations 
suffer from large uncertainties from a variety of sources: the large-x parton distributions of the 
projectile (x ~ 0.1 - 0.8), the large-z behavior of parton fragmentation functions (z > 0.6), NLO 
corrections to the "hybrid formalism" at high pj (see below), and so on. Presently, many eyes 
in this community are therefore focused on the upcoming p-i-Pb run at the LHC at yjs - 5 TeV 
which will open up phase space for small-x physics tremendously. This write-up shall therefore 
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also focus mainly (but not exclusively) on recent predictions for p+Pb collisions at the LHC 
based on the CGC approach. 

The frozen valence charge density per unit transverse area, p, in a hadron or nucleus is a 
random variable with a distribution determined by an effective action S [p\. In the limit of a large 
number of valence sources this distribution is Gaussian yj. 



r 2 

Suv[p\ = J ^ '''^y^ ■ 

The variance is proportional to the thickness of the target (= A^^^ for a nucleus, on average 
over impact parameters) since the sources act coherently. The classical field A^j is obtained 
from p by solving the Yang-Mills equations, analogous to the Weizsacker- Williams approach in 
electrodynamics. 

At the purely classical level the saturation momentum is energy independent. Energy de- 
pendence arises from quantum corrections to the classical field. It is not presently known how 
to incorporate all quantum corrections but a specific class of quantum fluctuations, those which 
are approximately boost invariant and are proportional to a, log 1 /x can be re-summed. For the 
two-point correlation function of light-like Wilson line operators, the so-called dipole scattering 
amplitude 

N(r)^^{t[V(0)VUr)) , (2) 
this is accomplished by the Balitsky-Kovchegov (BK) equation: 



dN{r,x ) 
d\n{xQl 



X) J 



ctn Kir, ri , rz) [TVCi , x) + N(r2, x) - N{r, x) - Nin , x) N(r2, x)] . (3) 



Due to the presence of non-linear effects this equation leads to saturation of the scattering ampli- 
tude at large dipole sizes r or small intrinsic transverse momenta kj- Evolution thus generates a 
dynamical scale Qs{x) which grows with energy (or 1 /x). In practice, eq. (O is nowadays solved 
with running coupling accuracy [5], and one of the main goals in this field is to test whether 
the evolution speed predicted by rcBK is consistent with observations or if improved accuracy is 
required. From the Fourier transform of N(r, x) (transformed to the adjoint representation) one 
obtains the (dipole) unintegrated gluon distribution (UGD) 0(k±, x) ~ k]_ NA(k±, x)las{k\). 

Solving (O requires an initial condition 7V(r, xq), where the reference rapidity log 1/xo is 
typically estimated to be Yq - logl/xo - log 100 - 4.6. N{r,xo) is the dipole scattering 
amplitude corresponding to the initial classical field at the reference rapidity Yq. In the MV 
model ([T]| it is given by 



N(r,XQ) = 1 - exp 



_i(,^e^(,„)fiog-l 



(MV: 7=1; AAMQS: 7 ^ 1.1), (4) 



with Q^(xo) ~ g^jj^ and 7=1. Ref. \^ performed detailed fits to (the most recent) HERA DIS 
data and found that a larger value for the initial "anomalous dimension" is preferred. Below, we 
shall see that semi-hard pj distributions in p-np collisions are consistent with the AAMQS initial 
condition coupled with rcBK evolution but clearly exclude the MV model initial condition for 
protons. As a first quantitative achievement we thus note that the proton rcBK-UGD at small 
X < xq (averaged over impact parameter) is now known to some degree of accuracy. However, 
the physical origin of the AAMQS correction to 7, which corresponds to a suppression of the 
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classical (^{kj^,x) ~ bremsstrahlung tail, is not very well understood. Since the valence 
charge density in a proton at - 0.01 is not very large one may consider corrections of order 
~ to the MV model action ij. 



2. Multiplicities and multiplicity distributions 

The p^-integrated multiplicity of charged particles is the most basic observable in particle 
collisions. For a variety of reasons though, the "details" are too many to discuss here, it is also 
hard to compute precisely, unfortunately. That said, one may hope that the main dependence 
of dN/drj on energy and system size is through the saturation scale Qs and so it is certainly 
warranted to check how this compares to data. 
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Figure 1: Charged particle rapidity distribution at 5 TeV from the KLN gluon saturation model 
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This issue was first addressed by the KLN model which indeed is in reasonable agreement, 
certainly well within its level of accuracy and credibility, with the centrality dependence of the 
multiplicity at midrapidity in Au+Au and d+Au collisions at RHIC, Pb-i-Pb collisions at LHC, 
and even p+p collisions from 0.9 - 7 TeV; for the most recent compilation and references to the 
original KLN papers we refer to ref. That paper also presented a prediction for p+Pb colli- 
sions at 4.4 TeV; an update for 5 TeV is shown in fig.[T] Other predictions for the multiplicity in 
p+Pb collisions at the LHC which instead use rcBK or IPsat UGDs, diff'er in their hadronization 
prescription or in the treatment of the nuclear geometry and of Glauber fluctuations, lead to very 
similar predictions at 77 = (deviations from the central value are + 15%) 

1 111 . This con- 
firms that the energy dependence of particle production is determined mainly by the growth of 
the saturation scale. It would be interesting also to compare theory and data for different cuts on 
A^part- The prediction is in very good agreement with preliminary ALICE data II12II at midrapidity, 
with a slightly too steep 77-dependence. 

One can also study the entire multiplicity distribution rather than just its mean. Assuming 
that particle production is due to strong classical fields, multiplicity fluctuations arise due to color 
charge density fluctuations of the effective sources at midrapidity. For a quadratic action such as 
in the MV model this leads (approximately) to a negative binomial distribution 111 311 : corrections 
are illustrated in ref. lll4ll . At midrapidity the bulk of the multiplicity distributions from pH-p 
collisions at energies 0.2 - 7 TeV can indeed be described quite well by NBDs [£, JJJ which 
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UA5 ppbar 200 GeV 
NBDfit<n>=3.9, k=1.5 
ALICE NSD pp 0.9 TeV 
NBDnt<n>=5.2, k=1.4 



ALICE NSD pp 2.36 TeV 
NBDfit<n>=6.8, k=1.2 
CMS NSD pp 7 TeV 
NBDfit<ii>=9.],k=l.l 
prediction MC-rcBK, p+Pb mb, Vs=4.4 TeV 
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Figure 2: KNO scaling plot of charged particle multiplicity distributions at \r]\ < 0.5 in NSD collisions at various 
energies |16] and NBD fits; z = N^\^I{N^\-^) and 'P(z) = (Wch) ^'(A'ch)- Note that the quoted mean multiplicity includes 
neutral particles; also, that here the fluctuation parameter k is integrated over the transverse plane of the collision. p+Pb 
prediction from ref |T3l . 



in fact exhibit KNO scaling to a good approximation lfl4l [Tsll . see fig. |2] Approximate KNO 
scaling has been predicted to persist even for minimum-bias p+Pb collisions at the LHC 015 1. 
This represents an important check of our current understanding of multi-particle production. 
Most importantly such intrinsic particle production fluctuations aff'ect the initial state for hy- 
drodynamics in A-hA collisions 1,15. XIl and thus could manifest in final-state flow and angular 
correlations. In the CGC approach, they are expected to occur on sub-nucleon distance scales on 
the order of ~ IjQs flT*], as shown in fig. [3] In this context it is interesting to note that CGC ini- 
tial state models which do not incorporate intrinsic particle production fluctuations |18] appear 
to be inconsistent with the distributions of angular flow harmonics presented at this conference 
by the ATLAS collaboration ||i9D. 
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Figure 3; Fluctuating energy density distribution at midrapidity and time t = 0.2 fm/c in an A+A collision [17 



3. Nuclear modification factor Rp\ipT) 



Transverse momentum distributions of produced particles provide more detailed information 
than pj-integrated multiplicities. In particular, they can probe the tail of the distribution of gluon 
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intrinsic kj. To compare to small-x QCD evolution one should restrict to pTl{{z)yls) e"^-' < xq ~ 
0.01 where (z) is the typical momentum fraction in fragmentation. 
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Figure 4: Left: Single-inclusive /7j--distribution in p+p collisions at 7 TeV [1 1.1: the curves show results obtained with 
rcBK evolved UGDs starting from two different initial conditions at xq = 0.01. Theory S'-factor fixed at p7-= 1 GeV. 
Right: nuclear modification factor SpPb(pr) at = (w/o rapidity shift) fllll . 

As a first step, it is important to check the spectrum in p+p collisions to establish consistency 
with the UGD fitted to HERA DIS data. Furthermore, a correct p+p limit is requked for trustwor- 
thy minimum-bias p+A spectra and for R^p, nuclear modification ratios. The result 111 ill is shown 
in fig. |4]for two different initial conditions. Clearly, the MV model initial condition with rcBK 
evolution is excluded by the LHC data which exhibits strong suppression of the high-fc^tails as 
compared to this UGD. On the other hand the AAMQS initial condition which has been carefully 
fitted to HERA data also agrees reasonably well with the CMS spectrum ||20|] and thus conforms 
to process independence of the (dipole) UGD. 

Fig.|4]also shows a prediction for the nuclear modification factor for p+Pb collisions at 5 TeV. 
The prediction and its many uncertainties are discussed in detail in ref. fTT). Generically, for all 
UGDs there is a suppression of particle production at = 1 GeV which grows stronger for 
more "central" collisions (higher A^part)- ^pPb then approaches and may even exceed unity at 
higher pj, as the evolution window shrinks; recall that in hadronic collisions, at fixed rapidity 
X oc pj. Also, fluctuations in the thickness of the target amplify higher-twist "anti-shadowing" 
effects. The shape of the Rppbipr) curves provides a test for the evolution speed predicted by 
rcBK; again, within uncertainties the prediction matches the preliminary data from ALICE. 

A very important recent development addresses corrections to the so-called "hybrid formal- 
ism" for particle production in asymmetric kinematic configurations (such as forward produc- 
tion in p+A collisions): large-x parton distributions are described within DGLAP while the 
small-x target field is obtained from rcBK, for example. Until recently only the LO expres- 
sion corresponding to elastic scattering of coUinear projectile partons on the dense target field 
was known (tA^. Ref. (l^ computed inelastic corrections corresponding to large-angle emission 
of a gluon in the projectile wave function which then exchanges little transverse momentum with 
the target. This correction is formally suppressed by one power of ff , but enhanced by log pr I Qs- 
Numerical evaluations 111 1112311 show that at high pj inelastic corrections drive /?pA up. However, 
for realistic values of as this contribution can be very large, especially for p+p collisions where 
Qs is relatively small 111]. The full NLO expression has been derived recently lf24ll and awaits 
numerical evaluation; this is crucial in order to understand the accuracy with which we are able 
to compute particle production in the forward region. 
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Figure 5: Nuclear modification factor RpPbiPr) at = 2, 4 obtained from the hybrid formalism llTIl 



A prediction for Rp^ at forward rapidities is shown in fig. |5] As expected, generically the 
suppression increases with rapidity and could grow much stronger than predicted by some current 
leading-twist shadowing approaches 1251,1261]. Notice the large correction (mainly to the p+p 
baseline) at high pr due to the inelastic term, even for very small = 0. 1 ; as already mentioned 
above more reliable predictions at high pr require fuU NLO accuracy. 

4. 7-hadron and hadron-hadron azimuthal correlations 

There has been a lot of activity recently also to predict azimuthal correlations which could be 
a very powerful tool for detecting high density effects. In essence, 2 — > 2 hard scattering should 
lead to approximately back to back correlations but if the projectile parton scatters off a strong 
target field that correlation should weaken or even disappear entirely. For the strongest effect, 
transverse momenta should be on the order of Qs, not far above. 

d + Au, 2.4 < 1/1,1/2 < 4, 1 GcV < p5f < pp'', p5?» > 2 GeV 



hadron: - 5 GeV. r\^= 3 

photon: k , - 3 
0.2 TeV __ k=3GeV 




A9 A<p [rad] 

Figure 6: Left: predicted y-hadron azimuthal correlations in the forward region of d+Au collisions at RHIC (27|| . Right: 
di-hadron coiTelations (incl. pedestal from double parton scattering) from a Gaussian approximation to JIMWLK evolu- 
tion (33l compared to (preliminary) STAR data [33j|. 



The most straightforward process for theory is y-hadron correlations. As only one colored 
particle interacts with the target, this process, within the hybrid formahsm, involves only the 
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DGLAP PDF of the projectile and the dipole unintegrated gluon distribution mentioned above 
(which can be obtained from rcBK). Theoretical uncertainties should therefore be reasonably 
small, allowing for quantitative tests of the theory. A first prediction for y-hadron correlations in 
the forward region of d+Au collisions at RHIC is shown in fig.|6][27]. 

Di-hadron correlations in the forward region are more intricate as now two colored particles 
interact with the target. Aside from the BK dipole the cross section for this process also involves 



at high pt, the cross section can be expressed in terms of the dipole and the (distinct) Weizsacker- 
Williams (WW) UGDs [22]. The evolution of the quadrupole with energy follows from the 
JIMWLK functional renormalization group equation. A numerical solution has been presented 
in ref. ['so'] and turned out to agree rather well, at least for a few simple configurations, with a 
Gaussian approximation (worked out in [29.] ) which expresses it as a (complicated) function of 
the dipole. A Gaussian approximation is extremely useful in practice since numerical solutions 
of the rcBK equation are rather straightforward; and because the dipole (incl. its initial condition) 
are rather well constrained by data, see above. The interest in correlations has also initiated new 
theoretical insight in form of a Gaussian approximation to JIMWLK evolution llJlil . 

To compare to data, ref. |32] also evaluates the pedestal due to double parton scattering. The 
authors observe broadening of the away side peak, fig.|6] due to the presence of non-linear effects. 
An earlier analysis using model parametrizations for the dipole and WW gluon distributions was 
presented in ref. ll34ll . 

5. Summary 

Tremendous progress has been made in recent years to develop the CGC effective theory into 
a quantitative framework which can confront data from RHIC and LHC. This includes (but is 
not limited to) improved computations of multiplicities for a variety of systems and energies; 
investigations of multiplicity distributions, KNO scaling in p+p, p+A and initial-state fluctua- 
tions in Ah-A collisions; semi-hard transverse momentum distributions and nuclear modification 
factors with rcBK UGDs and including Glauber fluctuations; photon production, y-hadron and 
di-hadron correlations. Currently, the CGC is the only formalism to address such a broad range 
of observables systematically and with some success. 

In the near future we should see further improvements of the accuracy and reliability of some 
of the predictions, for example a first computation of particle production in the forward region 
of p-Hp and p-i-A collisions at full NLO level. Moreover, there is great potential in understanding 
initial-state fluctuations in the "little bang", such as the scale on which they occur and how 
they reflect in the final state. Further, with respect to gluon correlations in the boosted wave 
functions which affect multi-hadron correlations. We are also presently lacking a description of 
the transition from semi-hard pj at small x to the DGLAP regime. These (and other) questions 
remain to be solved to achieve a more complete understanding of QCD, and the structure of 
matter, at short distances and high energies. 





(5) 
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